Effects of thickness on the spin susceptibility of the 2D electron gas 
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Using available quantum Monte Carlo predictions for a strictly 2D electron gas, we have estimated 
the spin susceptibility of electrons in actual devices taking into account the effect of the finite 
transverse thickness and finding a very good agreement with experiments. A weak disorder, as 
found in very clean devices and/or at densities not too low, just brings about a minor enhancement 
of the susceptibility. 
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Spin fluctuations are believed to play an important role 
in the two dimensional electron gas (2DEG), near the 
apparent metal-insulator transition observed at low tem- 
perature in clean devices, with lowering the density 0. 
Indeed, it has been found that the application of an in- 
plane magnetic field, which polarizes the electron spin, 
suppresses the metallic conductivity[l|,|2j. This, together 
with earlier suggestions that metallic behavior in 2D 
should be accompanied by a tendency toward a ferro- 
magnetic instability^, has recently prompted a num- 
ber of experimental investigations of the spin suscepti- 
bility X s of the 2DEG|1 Ull 0, S H ES [Hi, which is 
generally found to increase [l2j in an appreciable man- 
ner when decreasing the density n. Similar behavior is 
found also on the theoretical side for a strictly 2DEG, ac- 
cording to the most recent (and most accurate) quantum 
Monte Carlo (QMC) results[l3j|. However the suscepti- 
bilities measured in different devices differ among each 
other and, with one exceptional , d° not agree with the 
theory^j. Evidently, details of the devices play a role in 
determining the properties of the 2DEG and should be 
accounted for by the theory, as we shall show below. In 
particular in experiments the electron gas (EG) (i) has 
a finite transverse thickness, (ii) suffers scattering by a 
number of sources (scattering which in fact determines 
its mobility), and depending on the system (iii) occupies 
one or two degenerate valleys. In this Letter we address 
points (i) and (ii) for one-valley systems, exploiting the 
available QMC data. We find that taking into account 
the finite thickness of the specific device quantitatively 
brings into agreement theory and experiment and recon- 
ciles measurements on different systems, whereas details 
of scattering sources play a minor role. Furthermore, we 
offer some comments on the effect of (iii) valley degener- 
acy. 

At zero temperature and at given number density n, 
the state of the 2DEG can be specified by the spin po- 
larization £ = (nj — ni)/n. The spin susceptibility 
Xs = (9C/9B)b=o, which measures the ratio of the in- 
duced spin polarization to an in-plane applied weak mag- 



netic field B, is readily shown to be inversely propor- 
tional to the derivative (9 2 -E(£)/c?C 2 )c=o, involving the 
EG internal energy E(£). In fact minimization of the 
energy per particle E(() + (gbHsB/2 with respect to 
( yields the condition E'(() = —g^iisE /2 from which 
{d(/dB) B =o = -(g b fi B /2)/E"(0) immediately follows. 
An estimate of the spin susceptibility can be thus ob- 
tained from the knowledge of the internal energy E((). 

The effect of thickness on the 2DEG can be cast, in 
the simplest approximation, in terms of a device spe- 
cific form factor F(q) modifying in Fourier space the 
2D electron-electron interaction v(q) — 2-Ke 2 jeq into 
v(q) = v(a)F(a) |l4|. Rather than performing new sim- 
ulations for each device we have estimated the effects of 
thickness on E(£) and hence on Xs in a straightforward 
manner resorting to perturbation theory. In fact, to the 
lowest order in Av(r) = v(r) — v(r), one has for the en- 
ergy per particle E(Q = E 2D (C) + A(C), 



A(C) = - / drAv(r){g 2D ((;r)-l], 



(1) 

with £/2£>(C) an d Q?. r> (C l r) the known energy [l^ and 
pair correlation function L 15j, respectively, of the strictly 
2D electron gas. The accuracy of the energy estimates 
obtained in such a manner has been checked a posteri- 
ori per forming selected simulations with the interaction 
u(r)[lg|. We have computed the effect of thickness in two 
cases, for a GaAs HIGFET0 and for an AlAs quantum 
well (QW)E3. In the first case the form factor is Il7i 



££ + 3 9_ 1[1 + £i-3 

86 + 8fe 2n b l ' 



(2) 



with b 3 = 487rm(,e 2 n*/e?i 2 and n* = rid + j^n. Here, 
is the depletion charge density in the device^, mb = 
0.067m e the band electron mass and e = 12.9 the average 
background dielectric constant. For the AlAs QW, on the 
other hand, the form factor can be written as 19] 
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FIG. 1: Spin susceptibility in the 2DEG. The thick full curve 
is the QMC prediction for a strictly 2D svstem[l^|. Experi- 
mental results are given by the thick dashed curve for a GaAs 
HIGFETpl and by various s ymb ols, corresponding to differ- 
ent samples, for AlAs QW's[l(j. The thin dashed and full 
curves add the thickness effect onto the QMC prediction^! 
for GaAs and AlAs, respectively. Finally, the thick and thin 
dot-dashed curves provide the QMC prediction without and 
with inclusion of thickness for AlAs, as obtained from the 
polarization field Bp (see Eq. ^ in the text). The arrow 
indicates the location of the MIT transition in AlAs. 

with a = 45l the width of the well[l(Tj. Once F(q) 
is known, it is a simple matter to evaluate A(£), us- 
ing Eq. m, from which the enhancement Xs/xp °f the 
spin susceptibility Xs on its independent-particle or Pauli 
value xp is immediately obtained as Eq(0)/E"(0), with 
Eq(Q = Ep (1 + C 2 )/2 the energy per particle of non in- 
teracting Fermions in 2D and Ep the Fermi energy. Thus 

Xs/XP = E F /E"(0). (4) 

Our main findings |2Cj|| are summarized in Figure ^ 
which shows a number of calculations and measurements 
of Xs/xp- This quantity is plotted against the 2D cou- 
pling parameter r s = U/Ep = l/y / 7mas to get rid of 
uninteresting details of different materials which sim- 
ply determine the effective Bohr radius as = h 2 e/mbe 2 
through the dielectric constant e and band mass mi,; 
above, U = e 2 y / ™/e denotes a rough estimate of the 
potential energy per particle and Ep — h 2 im/mh- 

The QMC prediction for the strictly 2D case|ljj (thick 
solid line), is between 30% and 50% off the experimental 
values for GaAs (thick dotted line) . The key result of this 
paper is that this significant discrepancy is quantitatively 
explained as an effect of finite thickness, as clearly shown 
by the thin dotted line, obtained via Eqs. QEJ) and Q). 
This conclusion is further strengthened by our explicit es- 
timate of the effect of weak disorder, due to background 
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FIG. 2: Difference between the QMC energies of the unpolar- 
ized and polarized phases for the strictly 2DEG (solid line) 
and with the thickness of the GaAs HIGFET (dashed line) . 

doping in the GaAs HIGFET @, which turns out to be 
negligible (see below). We emphasize that the parame- 
ters entering the form factor of Eq. @ reflect our knowl- 
edge of the real sample [l^. and they are not adjusted to 
achieve a particular value of the spin susceptibility. 

In view of the interest for a possible ferromagnetic in- 
stability in low-density 2D electron systems, the question 
arises as to whether thickness, which noticeably changes 
the spin susceptibility, also alters the stability range of 
the polarized fluid (26 < r s < 35) predicted by QMC in 
the strictly 2D case|13j. Figure shows that this stabil- 
ity window is only slightly shrunk (by < 2 in r s ) and the 
ferromagnetic instability pushed at slightly higher r s . 

For the AlAs device the situation is somewhat dif- 
ferent. In order to engineer a one-valley 2DEG with 
isotropic mass, the electrons are confined in a very nar- 
row QW, thus reducing the importance of finite thick- 
ness effects, but also boosting the possible influence of 
well width fluctuations [l^. The spin susceptibility, mea- 
sured with either the tilted field or the polarization field 
methods (filled and empty symbols in Fig. ^ respec- 
tively), turns out to be fairly close to the strictly 2D 
QMC value. As expected, the tiny width of the QW does 
not affect this value significantly, as shown by the thick 
and the thin full curves in Fig. ^ Therefore the small 
discrepancy between the QMC prediction (with or with- 
out finite thickness) and the experimental result points 
to a possible role of disorder. 

Before discussing our estimates of the effect of disor- 
der, however, it is worth to briefly comment on the two 
different techniques used to measure Xs/xp- Within the 
Landau-Fermi liquid theory, the susceptibility enhance- 
ment can be expressed in terms of the quasiparticle pa- 
rameters g* and m* as Xs/xp = 9* m * I ' 9b"mb, where mt> 
and gt are the mass and g-factor entering the hamiltonian 
describing the interacting electrons, which for electrons in 
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a device coincide with the band mass and g-factor. One 
of the experimental techniques employed to estimate \s 
is the tilted field method of Fang and Stiles j^, which 
allows the determination of g*m* from the analysis of 
the minima in the Shubnikov-dc Haas oscillations. The 
experimental results in Fig. ^ for the HIGFET and part 
of those for the QW's (full symbols) were obtained with 
this technique. An alternative manner to extract Q*m* 
from experiments has been suggested by Okamoto|4j. If 
the interacting electrons can be replaced by independent 
particles with effective parameters (m* and g*), then the 
(in-plane) magnetic field necessary to induce full spin po- 
larization satisfies HBg*Bp = 2E F — h 2 2Tm/m* , which 
gives Xs/xp = ZEp/gbmbBp. However, the (in-plane) 
polarization field must also satisfy the exact condition 
gb^BBp = 2E'(l), which combined with the above yields 

Xs/xp = E F /E'{\). (5) 

The experimental results for the AlAs QW's obtained 
with both techniques are consistent with each other, due 
to the spread in the data. On the other hand, this is not 
the case with QMC 13] for which Eq. 10) yields an ap- 
preciable overestimate of the susceptibility enhancement. 
While the correct definition of the spin susceptibility is 
the one of Eq. Q , in comparing measured and calculated 
values of Xs/xp it is appropriate to refer to theoretical 
estimates consistent with the adopted experimental de- 
termination. In particular for low density, where only 
polarization field data are available (say r s larger than 
about 6), the theoretical value to be considered is the 
dash-dotted lines, calculated using Eq. (J5J. 

We now turn to the discussion of the disorder effects. 
A realistic description of these devices needs the inclu- 
sion of the different scattering sources that determine the 
mobility at zero temperature. The GaAs HIGFET 9] is a 
very clean device with no intentional doping and with a 
concentration of background impurities of ss 5T0 12 cm~ 3 , 
which is indeed a very low value [2 2j. Such a concen- 
tration has obtained through a best fit of the measured 
mobility [2^|, as a function of the electron density, com- 
puted using Born approximation}^ . We estimate the 
gross effect of such weak disorder on the spin suscepti- 
bility by means of perturbation theory, describing the 
disorder in terms of an external one-body potential u(r), 
coupling to the electron density, with known first and 
second moment ensemble averages |25|. As usual we as- 
sume a vanishing first moment. The first non vanishing 
contribution to the energy reads: 

with u(q) the 2D Fourier transform of the one-body po- 
tential and (• • ■ ) denoting the ensemble average over dis- 
order configurations, per unit area. Above, Xn,n{l,Q) is 
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FIG. 3: Effect of a weak disorder appropriate to the GaAs 
HIGHFET on the spin susceptibility. Upper and lower dotted 
curves are for the strictly 2DEG and for a 2DEG with the 
GaAs HIGHFET thickness. Other curves as in Fig. □ 

the density-density response function (in Fourier space) 
at polarization £ for the strictly 2DEG|24|. 

In Fig. El we show our results for the GaAs HIGFET, 
using Eq. JfjJ to estimate the effect of disorder in a density 
range in which A^(C) is much smaller (less than 1%) than 
the unperturbed 2DEG energy . The effect of disorder 
is to enhance the spin susceptibility, contrary to that 
of transverse thickness. However, in the experimental 
density range (3 < r s < 8) such an effect is negligible. 

The mobility in AlAs QW's is roughly two orders of 
magnitude smaller than in the GaAs HIGFET, as the 
scattering sources are more effective. We were able 
to reproduce the measured mobilities [Icj. using Born 
approximation !^ l2t| . only up to r s < 4. In this density 
range the additional enhancement of Xs due to disorder 
remains < 10%, yielding a good agreement between the 
prediction of Eq. 0. with £7(0 = £7 2 £>(C)+A(C)+A d (C), 
and the available results of the tilted field experiments. 
Our findings suggest that, far from the metal-insulator 
transition, occurring at r c s = 8.3 in AlAs QW's and at 
r c s = 12.4 in the GaAs HIGFET, disorder effects yield a 
small enhancement of Xs, which is either negligible as in 
GaAs or helps in reducing the small residual discrepancy 
between QMC and experiments. However, to obtain in- 
dications valid at larger r s and/or stronger disorder, an 
approach that takes into account disorder and electron- 
electron interaction on the same footing is required. 

We finally comment on the two-valley electron sys- 
tems, realized in Si MOSFETsji H H and in wide 
AlAs QW[lll. In particular Shkolnikov et al.ll], tun- 
ing the valley population, have shown that valley de- 
generacy brings about a depression of the spin suscep- 
tibility, in sharp qualitative contrast with the enhance- 



4 



ment predicted by Hartree-Fock theory. A first esti- 
mate of the spin susceptibility of a two-valley symmetric 
EG can be simply obtained from previous QMC studies 
of the energy of four-|27) and two-component [2^ elec- 
trons performed with the same level of accuracy, assum- 
ing a quadratic dispersion of the energy E(() with cH^I- 
Such an estimate clearly shows the qualitative effect ob- 
served in Ref. [ill A detailed comparison with either 
the Si MOSFEtUHHQ or the anisotropic-mass AlAs 
QW[ll| devices along the lines of the present calculation 
would require QMC input which is presently not avail- 
able. We are currently performing extensive simulations 
of the strictly 2D two-valley system with finite polariza- 
tion to obtain an accurate theoretical prediction of x s in 
such a system. 

In conclusion we have shown that a realistic descrip- 
tion of actual devices, starting from the 2DEG model and 
including specific features of the systems, enables us to 
reproduce the spin susceptibility without any adjustable 
parameters. In GaAs HIGFET 9] the thickness plays a 
crucial role, while the weak disorder provides only a neg- 
ligibly small enhancement. In the AlAs QW's case[10j 
the 2DEG spin susceptibility is directly comparable with 
experiments. A residual small discrepancy is likely due 
to the influence of QW's width fluctuations. 
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